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Abstract 

We investigate the conditions under which particle multiphcities in high energy 
colUsions are Boltzmann distributed, as is the case for hadron production in e~^e~ , pp, 
pp and heavy ion colhsions. We show that the apparent temperature governing this 
distribution does not necessarily imply equilibrium (thermal or chemical) in the usual 
sense, as we explain. We discuss an explicit example using tree level amplitudes for 
N photon production in which a Boltzmann-like distribution is obtained without any 
equilibration. We argue that the failure of statistical techniques based on free particle 
ensembles may provide a signal for collective phenomena (such as large shifts in masses 
and widths of resonances) related to the QCD phase transition. 
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1 Introduction 



The use of thermal or statistical models to describe multiparticle production has a long 
history [1]. Recent studies [2-12] observe that the multiphcities of hadron production in a 
variety of contexts (e'^e", pp, pp and heavy ion collisions) are extremely well described by 
models involving thermal distributions of free hadrons. The only free parameters in their 
analysis are the temperature and volume of the model thermal system, and a parameter 
reflecting the level of equilibration of strange particles (in the heavy ion case the baryon 
chemical potential is an additional parameter). One might conclude from their results that 
thermal and chemical equilibrium among hadrons is already reached in individual jets pro- 
duced in high-energy scattering. 

However, this conclusion is somewhat unjustified, given that any mechanism for produc- 
ing hadrons which evenly populates the free particle phase space will mimic a microcanonical 
ensemble, and therefore yield apparently thermal results^. It is important to remark that 
this type of apparent thermalization will always yield Boltzmann weights which are func- 
tions of the free particle energy, and hence describe a non-interacting ensemble. True QCD 
thermalization, of the type associated with the quark-gluon plasma (QGP), and that exper- 
imentahsts hope to observe at RHIC, involves large collective effects, and hence is probably 
poorly modelled by a non-interacting hadron gas. 

Here we address the problem of deducing whether a process which leads to multiparticle 
production is thermal. We make the important distinction between phase space dominated 
phenomena, which lead to ensembles governed by free particle Boltzmann weights, and 
interacting thermal ensembles, in which collective effects can be important. We argue that 
a process which generates data that can be fit using free particle ensembles is merely good 
at populating phase space in a uniform way - it has not necessarily produced an interacting 
thermal region. In QCD, a region of this type with temperature of order 100 MeV or 
more should exhibit strong collective phenomena, such as hadron mass shifts, that preclude 
description in terms of a free particle ensemble. Hence, we argue that the failure of statistical 
techniques based on free particle ensembles should be regarded as a signal for the onset of 
true equihbration at heavy ion coUiders such as RHIC. 

This paper is organized as follows. In section 2 we review the relation between micro- 
canonical and canonical ensembles in statistical mechanics. We argue that multiparticle 

^After completion of this work we learned that a similar conclusion was reached by C. N. Yang et al. 
in [13]. They refer to the apparent temperature as a partition temperature and stress that no thermal 
equilibrium is implied. We thank Professor K.S. Lee of Chonnam National University in Korea for making 
us aware of this earlier work. 
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production in many cases is equivalent to a method of populating a modified microcanonical 
ensemble. Such an ensemble will produce "thermal" behavior even if there is no subsequent 
interaction of particles once produced. In section 3 we apply our results to multiparticle 
production and examine a specific example in a toy model involving tree-level photon pro- 
duction. In this toy model there is clearly no real thermal equilibrium - once produced, 
the photons do not interact. However, a quasi-Boltzmann distribution results in the limit 
where a large number of photons is produced. In section 4 we discuss the implications of 
our conclusions for heavy ion collisions. 



2 Microcanonical vs Canonical Ensembles 

Here we give a brief review of the relationship between the microcanonical and canonical 
ensembles (MCE and CE, respectively) in statistical mechanics. Recall that the MCE sums 
only over states with some fixed total energy, while the CE sums over all states with Boltz- 
mann weight e~^^ . The main result is rather familiar: under certain general assumptions, 
quantities computed in the MCE differ from those computed in the CE by an amount which 
vanishes as the size of the system N is taken to infinity. The importance of this result is as 
follows: the cross section for particle production in a high energy collision can be written as 
the expectation of the matrix element squared in the MCE corresponding to the free theory. 
In the limit of large N, one can therefore rewrite the usual phase space integral appearing 
in a cross section in terms of the average of the matrix element squared in a CE, which is 
controlled by the Boltzmann factor. This naturally leads to a certain "thermal" behavior 
which we discuss below. 

In the MCE the total energy is fixed and the probability density is constant over all of 
phase space. The computation of the entropy S{E) is as follows: 

r{E)=J26{E,-E), (1) 

s 

S{E)^\nr{E). (2) 

In Eq. (1), Eg is the energy of the state s. It is often more convenient to replace the delta 
function in (1) with the factor A{E) = exp{—(3{Es — E)), yielding a new quantity, the CE: 

f(£;) = ^e-^(^^-^). (3) 

s 

In the canonical ensemble there is no restriction on the energies of the states s. However, 
they appear with Boltzmann weight ex.p{—(3Es). The new quantity introduced, /3, describes 
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the temperature of the system, which is fine-tuned to ensure that the average energy is E. 
To see this, rewrite T{E) as follows: 

= J dE' T{E') e-f"^^'-^^ 

= y e^(^')-^(^'-^). (4) 

Now evaluate (4) in the saddle-point approximation. (This is also known as the Darwin- 
Fowler method [13].) Let 



dE' 



(5) 



so 



f ^ g5(£) J ^^f ^lS"{E){E'-Ef+- 



V ~S"{E) e + ■ 

Here we have assumed that S'{E) > (positive temperature) and S"{E) < (positive 
specific heat). It is easy to see that the difference between the CE entropy S{E) = lnr(£') 
and the MCE entropy S{E) is of order (In N)/N. The terms represented by the ellipsis in 
(6) lead to even smaller corrections, and will be neglected. 

Now consider a generic operator O. The average in the CE is given by 



e 



siE)^^l^ J dE' e'm-P{E'-E) ^o)m{E'), (7) 



where {0)m is the average taken in the microcanonical ensemble. Its logarithm can be 
expanded as follows 

H{0)m{E')] = H{0)m{E)] + (E'-E) ]n'[{0)M{E)] + ^iE'-Ef ]n"[{0)M{E)] + • • • . (8) 

The integral in (7) can again be performed in the saddle-point approximation. Note that 
if the operator O is of order N (as in the case of a particle multiplicity), the coefficients 
appearing in the expansion (8) are of order In A^. They lead to small shifts (of order (In N)/N) 
in the saddle-point value of the temperature /? and the overall prefactor. Thus the canonical 
and microcanonical averages of particle multiplicities converge in the limit of large A^. 

In the next section we apply these results to the computation of cross sections for particle 
production in high energy collisions. 
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3 Multiplicity Results 



The cross section for the process A + B ^ n particles is 

First, let us assume that the function^ |7Vlp/n2-E/ in the dominant region of phase space 
is slowly varying (we will relax this assumption shortly). Then, we can use the following 
approximation 

1 \M\^ 



2Ea2Eb\va-vb\ Uf'^E 

where M., E are averaged quantities. As noted, the integral in (10) is just the microcanonical 
ensemble for N free particles, and hence leads to "thermal" properties of the particle distri- 
butions and multiplicities. Considering the more general case where the number of particles 
is not fixed, we simply sum over all cross sections, 

a^^aAB^n, (11) 

n 

to obtain a MCE without fixed particle number. In the usual thermodynamic limit this sum 
is dominated by some particular value of n, so it is equivalent to consider the earlier case 
with n — N. 

We now treat the matrix element more carefully, by retaining it in the phase space 
integral. The resulting integral can still be turned into a canonical ensemble using the 
result of the previous section, provided the modified "entropy" {i.e., the logarithm of the 
phase space integral including the matrix element) continues to grow with total energy E 
(positive temperature), and the second derivative of this entropy with respect to E is negative 
(positive specific heat^). The energy-momentum delta function can then be replaced by a 
Boltzmann weight exp[(p^+j9B— X]p/) ■ which in the center of mass frame reduces to 
exp[/?(£^— Y,Ef)]. This yields the following result for the differential cross section: 

(fa 



dpi 



cx 



i- e-^'/- / n d^, -1 \M\' exp (- E E,/t\ . (12) 



^In relativistic field theory, we usually adopt the normalization convention {p\p') = 2Ep {2'jt)^6^{p —p'), 
which leads to the factor of l/2Ep in the phase space density. On the other hand, in statistical physics, 

wc usually adopt an cncrgy-indcpcndcnt normalization. Although the final result for the cross S(K;tion is 
independent of our normalization convention, what wc mean by "phase space dominated" depends on what 
we choose to be the unit of phase space. In this paper we will always be referring to the statistical mechanical 
unit of phase space, i.e. d^p. 

^These requirements are satisfied in the toy model we consider below, where the modified entropy behaves 
as S{E) ~ Ai'lnln(£;/m), for N "photons" with total energy E. 
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Note the natural appearance of the Boltzmann factor in (12). 

To proceed further, we need a specific model for the behavior of the matrix element. 
There are very few cases in which the matrix element for iV-particle production is explicitly 
known. One such case is the QED process 

// — > 1 spin up photon + {n — 1) spin down photons, (13) 

where the (massless) fermions have opposite spin. The matrix element squared for (13) is [14] 

where P, Q and Pi are the momenta of the two incoming fermions and n outgoing gauge 
particles respectively. In the (P, Q) center of momentum frame this may be written as 

\MiP^, Qi; It, 2;, • • • , nO\' = ]CEl{l + cos 0,)' f[ \ . , (15) 

i=2 Sill 

where E^, 9i are the energy and production angle of the ith photon, and we have lumped all 
of the remaining constant factors into AT. Using this to solve for the differential cross section 
yields {i ^ 1) 

The resulting number density of (spin down) photons is then 

""^/Ir^ra*'"'' '''' 

In this particular example we have a problem, because the photons are massless and there 
is an infrared catastrophe due to arbitrarily soft photons. Of course, the number density 
of observable photons - those with some minimum energy and angular separation from 
the initial // pair - is finite. For our purposes, we can always eliminate this problem by 
introducing a photon mass by hand in our toy model. In fact, we can introduce several 
species of "photons" with masses mi. Then, the abundance of each species is given by^ 

y (27r)3 (p2 + ^2)3/2 e ■ U»j 

The integral in (18) differs from the one appearing in a pure Boltzmann distribution due 

to the factor of (p^ + mf)~^/^. Without this extra factor, the integral simply reduces to 

mfTK2{mi/T), where K2{x) denotes the modified Bessel function of order 2. 

^In writing Eq. (18) we have chosen to ignore the angular dependence. Had we chosen to retain it, the 
angular integration would produce the factor ln[(l + ^/^ — 4mf/s)/{l — -\/l — Arnj/s)] , where ^/s is the 
center of mass energy of the collision. 
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Figure 1: Multiplicity vs. particle species mass of the toy-model, and thermal fit. 



In a sense, the additional factor makes only a small difference relative to the exponential: 
when taken into the exponent it is of order hi{f3E), versus PE for the Boltzmann factor. 
However, actual ratios of particle abundances fii/rij will differ from thermal ratios. In Fig. 1 
we show the result of the multiplicities from (18) and a thermal best fit. While we didn't 
use any additional parameters, such as individual chemical potentials, the eventual quality 
of the fit would probably not be as good as what is observed in e'^e~, pp, pp and heavy ion 
collisions [2-12]. In other words, the hadronization process probably populates free particle 
phase space somewhat more evenly than our toy model. However, our toy model does 
demonstrate that Boltzmann-like distributions are not necessarily indicative of real thermal 
(chemical or kinetic) equilibrium. 



4 Discussion 



In the previous sections we argued that multiparticle production can readily lead to thermal 
behavior if the the process in question is phase space dominated. Because phase space 
is determined by free particle kinematics, the results correspond to an ensemble of non- 
interacting particles. In other words, because the arguments of the energy delta function 
in the cross section Eq. (9) are simply free particle energies, the corresponding Hamiltonian 



7 



appearing in the ensemble in Eqs. (1) and (3) is the free Hamiltonian, with no interactions. 
Our toy model of photons suggests that this result is rather generic in any process where a 
large number of particles is produced. 

If the ensemble is dominated by a particular species of particle {i.e. the lightest particle), 
the apparent "temperature" will be related to the mass m of that species. This is because 
phase space is maximized by producing as many particles as possible, each with a kinetic 
energy of order m. Suspiciously, the typical temperatures produced by the excellent thermal 
fits of e"'"e~, pp, pp and hescvy ion collisions [2-12] are all of order the pion mass. 

Of course, an ensemble of non-interacting hadrons is not very interesting. It gives us 
no information about the actual QCD phase diagram. In real QCD the energy of a state 
consisting of many particles is modified due to interactions: it is not simply the sum of the 
free particle energies. At high density or temperature interaction effects are large and lead 
to large deviations from free particle results. Were this not so one could never see collective 
phenomena such as chiral symmetry restoration or a deconfinement phase transition. In an 
ensemble of interacting hadrons we already expect significant collective effects at tempera- 
tures T ~ m,r, such as a decrease in the value of the quark condensate {qq)T, and shifts in 
the various hadron masses. 

We expect these effects to lead to the failure of statistical techniques based on free particle 
ensembles, which predict that the multiplicities should fall roughly exponentially (as e^"'/"^, 
see Fig. 1). However, if the masses and widths are shifted from their vacuum values at 
the instant of chemical freeze-out, as would be expected at the ~ 170 MeV temperatures 
obtained in the fits, then a plot of multiplicity versus vacuum mass will deviate from the 
thermal prediction: i.e. it will have the form^ q-'^{t)/t ^ ^j^g other hand, it is also possible 
that the system will remain in equilibrium long enough for the masses and widths to return to 
their vacuum values. If so, then a thermal fit will perform well, but the resulting temperature 
would not be near the chiral phase transition. 

Knowledge of these mass shifts would be necessary for any detailed predictions. Although 

there have been some promising results [15] . it has proven difficult to extract this information 

from lattice data. A perturbative estimate of the decrease in the quark condensate was 

derived in Ref. [16]. In the chiral limit (m„ = = 0), an analytic result is possible. It 

■''Wc have simplified our discussion here in two ways. First, we have ignored the possibHty of introducing 
chemical potentials related to conserved quantities. Since the number of data points to bo fit is much larger 
than the handful of these potentials which may justifiably be introduced, it is highly unlikely that the effects 
of all of the mass shifts could be reproduced in this manner. Second, many of the hadronic states decay 
rapidly and so affect the relative populations of the observed particles. Inclusion of these effects adds no 
degrees of freedom to the fits. It would be amazing if the shifts in masses and widths should conspire to 
reproduce the vacuum results. 
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reads 

{m)T = (gg)T^o [1 - ^ - § (l^ J + ^(^')J ' (^9) 

where F ^ 93 MeV is the (zero temperature) pion decay constant. Using real world quark 
masses changes this result only slightly. The corresponding estimate for the temperature 
of the chiral phase transition {i.e. the point at which the condensate is essentially zero) is 
170 MeV in the chiral limit, and 190 MeV in the real world, consistent with lattice data. 

We expect the thermal pion mass to obey the finite temperature version of the Dashen 
formula 

ml{T) c -^|(g-g)T|, (20) 
where (again, for two massless flavors) [16] 

r2 



F\T) = F\0) 



6F2(0) 



+ O(T^) 



(21) 



Note that at leading order the thermal pion mass actually increases slightly with temperature. 
The temperature dependence of the baryon masses is a harder problem, but in the naive 
quark model we might expect them to behave roughly as 

mB{T) ~ 3 \{qq)T\, (22) 

where the constituent quark mass is simply due to the condensate. This estimate at least 
incorparates the fact that the baryons must become nearly massless at the chiral phase 
boundary, since chiral symmetry prevents them from obtaining a mass. Equations (20) 
and (22) suggest that the relative abundance of baryons compared to pions will increase at 
high temperature. While we don't necessarily beheve that (22) is very accurate, the point is 
that the thermal pion and baryon masses probably do not depend on the quark condensate 
(and hence the temperature) in the same way. Thus, it is probably inconsistent to imagine 
fitting the properties of an interacting hadron gas at T ~ using vacuum hadron masses. 
Yet, essentially all recent heavy ion data^ agrees well with multiphcities generated by free 
thermal models, with temperatures of roughly T ~ 50 — 170 MeV. 

In Fig. 2 we reproduce a plot from the paper of Cleymans and Redlich [12], where the 

fitted temperatures and chemical potentials resulting from LEP, CERN/SPS, BNL/AGS and 

GSI/SIS data are displayed. One interpretation of these results (which we do not subscribe 

to) is that a large region of the QCD phase diagram in the temperature-density plane has 

already been explored! There is no question that the quahty of the free thermal fits is 

^In fact, from this point of view any model (such as a parton cascade model) which reproduces free 
particle thermal multiplicites at temperatures of order 150 MeV probably lacks some important dynamics 
associated with the phase transition. We might classify it as just another efficient populator of phase space. 
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Figure 2: Temperatures and chemical potentials obtained by Cleymans and Redlich [12] 
from thermal multiphcity fits of LEP, CERN/SPS, BNL/AGS and GSI/SIS data. 

quite good. However, this in itself suggests that an interacting thermal region has yet to be 
produced in these experiments. Rather, it is quite possible that the collisions simply serve 
as a mechanism for populating phase space, without ever evolving through configurations in 
real thermal and chemical equilibrium {i.e. actual points on the phase diagram in Fig. 2). 
If the system had passed through real equilibrium, we suggest that the observed final state 
multiplicities could deviate significantly from those which can be generated by free thermal 
models. Thus the failure of such models to fit the data could be a signal for real equilibrium. 
This idea is explored in the context of phenomenologically motivated models in [17]. For more 
recent work related to this paper, which originally appeared as preprint nucl-th/0001044 in 
2000, see [18]. 
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